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ICASE Fluid Mechanics 


Due to increasing research being conducted at ICASE in the field of fluid mechanics, 
future ICASE reports in this area of research will be printed with a green cover. Applied 
and numerical mathematics reports will have the familiar blue cover, while computer science 
reports will have yellow covers. In all other aspects the reports will remain the same; in 
particular, they will continue to be submitted to the appropriate journals or conferences for 
formal publication. 
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and 

Demetrios T. Papageorgiou 1 
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Newark, New Jersey 07102 


ABSTRACT 

We investigate the stability of two phase Couette flow of different liquids bounded between 
plane parallel plates. One of the plates has a time dependent velocity in its own plane, which 
is composed of a constant steady part and a time harmonic component. In the absence of time 
harmonic modulations the flow can be unstable to an interfacial instability if the viscosities 
are different and the more viscous fluid occupies the thinner of the two layers. Using Floquet 
theory, we show analytically in the limit of long waves, that time periodic modulations in 
the basic flow can have a significant influence on flow stability. In particular, flows which 
are otherwise unstable for extensive ranges of viscosity ratios, can be stabilized completely 
by the inclusion of background modulations, a finding that can have useful consequences in 
many practical applications. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-19480 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (1CASE), NASA Langley Research Center, Hampton, VA 23681. 
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1 Introduction 


It is well known that plane Couette flow of two superposed fluids of different viscosities 
can be unstable. This instability is absent when the viscosities are equal and is a result of 
an interfacial deflection growing due to viscosity stratification. This mechanism was first 
described by Yih (1967), who used perturbation theory to obtain an analytic expression for 
the growth rate in the long wavelength limit. The results can be summarized as follows (for 
clarity, in the absence of gravity): Instability is possible only if the viscosities are different. 
Arrangements where the more viscous fluid occupies a thinner layer than the less viscous one 
are unstable, whereas converse arrangements are linearly stable to long waves. The instability 
is present at all Reynolds numbers, although in the limit of zero Reynolds number the growth 
rates become asymptotically small. Flow, therefore, is necessary to excite growing waves. 
Moderate surface tension is negligible to leading order for long waves, but acts to stabilize 
very short waves. Yih’s work was extended more recently by other investigators. Hooper 
and Boyd (1983) considered the linear stability of two co-flowing viscous liquids of different 
viscosities separated by an interface and extending to infinity. This is Yih’s problem in the 
absence of walls and can be used as a useful model in the study of short waves (short waves 
at the interface would not feel the effect of the walls, to leading order at least). Instability 
is found in all cases. Long waves are stable for all viscosity ratios, while short waves, in 
the absence of surface tension are unstable with asymptotically small growth rates - when 
included, surface tension damps out all waves which are short enough. Disturbances with 
general wavelengths are unstable and a maximum growth rate is attained at some viscosity 
ratio. This work was extended by Hooper (1985) who included a bounding moving wall for 
the lower fluid. The flow is shown to be unstable to long waves as long as the bounded layer 
is also the more viscous one, in agreement with Yih’s findings. An explicit expression for 
the eigenvalues is also given in terms of Airy functions. 

The stability of Yih’s problem for general disturbances and including density differences 
and surface tension, was first solved numerically by Renardy (1985). She presents specific 
examples that show that the general long wavelength result described by Yih and Hooper 
and Boyd, can be extended to interfacial waves of arbitrary wavelengths. Another interesting 
feature is the finding of a second mode of instability at relatively high Reynolds numbers; 
this second mode has long wavelength and growth rates comparable to the first interfacial 
mode. Renardy (1987) carried out an analytical study of the same problem but for fluids 
with slightly different mechanical properties. Analytical expressions for the eigenvalues are 
constructed and instability is established in the thin layer limit if the thin layer is occupied 
by the more viscous fluid. 
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The stability characteristics of oscillatory flows are much more complicated than those 
of their steady counterparts. There are many examples where the inclusion of an oscillatory 
component to the steady flow can enhance or reduce stability (for plane Poiseuille flow, for 
example, see Grosch and Salwen (1968), Hall (1975), von Kerczek (1982); for a review on the 
stability of oscillatory flows see Davis (1975)). As far as we know the stability of oscillatory 
two-phase viscous flow has not been studied. Yih (1968) considered the stability of a viscous 
fluid layer on a flat plate performing a simple harmonic, motion in its own plane. There is no 
viscosity stratification (the upper fluid is air) and the flow in the absence of the oscillation 
contains no steady velocity component and is linearly stable. Yih showed, using a long wave 
Floquet theory, that the oscillatory flow can become unstable even though in the absence 
of oscillations the flow is stable. This was extended by von Kerczek (1987) to flow down 
a vertical plate which is performing a simple harmonic motion in its own plane. This flow 
is unstable even in the absence of oscillations and von Kerczek uses Yih’s long wavelength 
expansion to establish windows of instability. Our interest is to apply such a study to an 
oscillatory two-phase flow. In the problem we study here, depending on the fluid occupation 
areas in the unperturbed state (see earlier comments), the flow can be stable or unstable. In 
the former case, therefore, a quasi-static approach would yield background velocity profiles 
for each instant in time (time is treated as a parameter in quasi-steady theories) which are 
linearly stable while in the latter instance instability ensues for all parametrized profiles. A 
more yielding approach is that of Floquet theory where stability or instability is judged on 
the overall growth or decay with time of a perturbation over a complete period of the forced 
oscillation. Analytically this means that time-periodic solutions are constructed which can 
amplify, remain neutral or decay depending on whether the Floquet exponent is positive, 
zero or negative. We present representative results for several cases which in the absence 
of background oscillations the flow is unstable while imposed oscillations can stabilize the 
flow completely. On the other hand, flows which are stable can become unstable. We note 
that all the results given here are in terms of explicit but long formulae which yield stability 
results with little computational effort. Such explicit solutions are of considerable value in 
the generalization of the stability problem to arbitrary wavenumbers by use of continuation 
methods for instance. 

The article is organized as follows. Section 2 derives the undisturbed flow as an exact 
shear flow solution of the Navier-Stokes equations. In Section 3 the linear stability problem is 
formulated and the interfacial conditions are written down. Explicit solutions are developed 
by an expansion procedure carried out to three orders, so that we may determine the first 
stage where the expansion of the Floquet exponent is a real term. In Section 4 we present 
the results of the stability analysis and finally in Section 5 we draw some conclusions. 
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2 The basic state 


Two incompressible fluids of equal densities, p, and different viscosities, occupy a region 
of depth L between parallel infinite flat plates. The lower plate is fixed and the upper 
plate moves in its plane, with a steady velocity Uq , together with superposed sinusoidal 
oscillations, so that its velocity at time t* is given by C/ 0 + A cos (u>t*). The fluids are assumed 
to be immiscible and form separate layers, the upper fluid has viscosity p. i, the lower fluid 
has viscosity p -2 and a depth D. We denote these as regions (I) and (77) respectively. 
We look for exact solutions of the Navier-Stokes equations for a constant interface position 
and velocity component in the streamwise direction alone which depends on time and the 
vertical coordinate only. Using cartesian coordinates (x*,y*), the exact flow is described by 
the following equations 


dF[ ' m d 2 u\ 
dt* ~ P d y * 2 ' 
dUl = p 2 d 2 JT 2 

dt * P dr 2 ’ 


L>y*>D 
D>y*> 0 


The boundary conditions are, no slip at the walls, 


(0, 

(//). 


U x {y‘ = L) = Co + A cos (uit’) , 

P* 2 ( y* = 0) = 0. 


and continuity of velocities and stresses at the interface, 


TTAy-^D) = Z^(/ = C), 

"'ly (s ' = i,) = ^ (y ' = D) - 

The solution is easily found by separation into a steady part and an oscillatory part of the 
same frequency as the plate motion. The basic flow is, 

jj . = wlV - C.Pfe - £t) + AR (U tav) + e i-'"'-)}, 

f.iL* — D ( #/« — flil v- L j - J 

'rr) , 


tr 2 = 


Lp 2 - D(p 2- di) 

P\U 0 y 


Lpi - D(p 2 - /ii) 


+ A.7Z < 2K sinh 


m 2 


where IZ denotes the real part, ft* = (ipu>/pi)* and we define the ratio of the lower fluid 
viscosity to that of the upper layer by the parameter m = p-ij p \ . The constants K , L\ and 
L -2 are given by 

1 


I< = 


2 sinh cosh(/3* ( D — L )) + 2 m* cosh ^ sinh (ft* (D— L)) 


( 1 ) 


3 



( 2 ) 


L\ 

h 


K exp (—(3*D) 

r 

K exp (/?* D) 


in 2 cosh j 

!3*D 


Wsinkf^ 


sinh 


m 2 / 

i — m 2 cosh 


m 2 


\ m 

i - 

l. ma 


?)]• 


( 3 ) 


The flow is non-dimensionalized next with D as lengthscale and Uq as velocity scale, while 
time t* and pressure P* are non-dimensionalized by D/Uo and pUl respectively. This intro- 
duces the following non-dimensional groups (used later also): 


a = A ? = itd, 

I'O 

oR\ _ uD 

pDUg' “ U 0 ' 

The parameter a measures the relative depths of the two fluids, A measures the magnitude 
of the externally imposed oscillations, (3 is a measure of the Stokes layer thickness due to 
the oscillatory flow, R e is the Reynolds number, and fl is the non-dimensional frequency of 
the plate oscillations. Under these scalings the non-dimensional no-slip condition becomes 


a = — 


R, = 


L 

D' 

U 0 Dp 


U\ (y = a) = 1 -f A cos (Of) , 

U 2 (y = 0 ) = 0 , 

and the basic flow is: 

_ = my-m+1 +An ( \ Llt m + L 2 e^ 0y) ] e (,nt) } , 

am — 771 + 1 L 1 J J 

U 2 = + AT l {‘2K sinh ( e ( ‘ ni) l . 

am - m + 1 ( \m 2 / J 

The constants K, L x and L 2 appearing above are in their non-dimensional forms which are 
readily available from (l)-(3) by replacing (3*D by (3. When A = 0 we recover the basic 
Couette flow of different fluids in a channel (see Yih (1967)). 


3 Stability equations. 

The equation governing the linear stability of parallel shear flows is the Orr-Sommerfeld 
equation and arises from linearization of the Navier-Stokes equations and a normal mode 
expansion. In the present problem the background flow is time-dependent and so the eigen- 
functions depend both on t and y. In particular, the perturbation streamfunction is taken to 
have the form + 2 (x, y , t) = <f> i, 2 (y, <) e iax where a is the wavenumber of the disturbance and 
subscripts 1,2 denote quantities in regions (I) and (II) respectively. The stability equations 
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in regions (I) and (II) are, then, 


'a ^va 2 

Ft +mU ', 



A , . , a 2 (/, . 

a I <t>\— ia«t>i 


— a 2 ] (f>2 — ioc<f>2 


W 

d 2 U 2 

dy 2 


Re 


^ 4 01 9 2^01 , 4 / 

-2a‘ s ^rrr + a 0i 


m 

& 


dy 4 

d A h 


dy 4 


dy 2 

2 + ^02 I 


-2«^ 


dy 


(4a) 

(4b) 


The linearized boundary conditions are 

■Ma) = 0=^(0), (5) 

h(0) = 0 = ^ (0) . (6) 

We now let the interface position be y = 1 -f Sh (t)e' ax where S « 1 is the infinitesimally 
small amplitude of the perturbations. The linearized interfacial conditions are found by 
expansion of the exact conditions about y = 1 in powers of 6 and retention of the leading 
order contributions: 


0 i = 02, 

00i 002 


d 2 <j) 


dy 


dy 

r + <* 2 <f > i 


dy 


(1 - to) dUi 
m dy ’ 


m 


+ o: 2 02 


** r _Sa»*L_2«>*. W ' 


= m ^ 


dy 3 

d 3 0 2 

dy 3 


dy 


dy 


- Re 




dy 


dy 


- Re 


+ h 

L d y 2 

d 2 0i , . tt 00i • , 

, ac 2 ' 


+ 


ic^h’j 
Re ‘ 


The linearized kinematic condition at y = 1 reads 

dh , , _ _ , , » 

— + iahU\ + ia<pi = 0 . 

at 


(7) 

( 8 ) 

(9) 


( 10 ) 


( 11 ) 


Equations (4-11) constitute the partial differential system together with boundary condi- 
tions that governs the stability of the flow. The problem contains time and one space variable 
and in general requires integration in time as an initial boundary value problem. Stability or 
instability is determined by the large time evolution of initial perturbations corresponding 
to different wavenumbers. 

In this work we proceed analytically by studying the stability characteristics of long 
waves (a << 1). Such analysis provides a significant amount of information since it enables 



development of closed form expressions which can be used to check whether unsteadiness in 
the background flow can affect the evolution of linear perturbations. In the Couette problem 
studied by Yih (1967), instability was established for all Reynolds numbers by an analogous 
long wave expansion. The analysis of long waves is also useful in guiding numerical solutions 
of the initial value problem at general wavenumbers. 


3.1 Solution for long waves 

We now proceed analytically by considering the asymptotic limit a — * 0 which corresponds 
to disturbances with long streamwise wavelengths. We use a method similar to the Floquet 
theory applied by Yih (1968) to obtain analytically the growth-rate of the disturbance at 
the interface. We look for solutions for the streamfunction and interfacial amplitude as a 
power series in a <C 1 of the form, for j = 1,2, 

4>j ( y , t) = {<j>jo ( y , t ) + a<j)ji ( y,t ) + . . .) exp ((#o + a &i + • • -)0 ex P (ikott ) , 
h(t) = (h 0 (t) + a/ii (<) + •••) exp ({9 0 + a9 x + . . .)t) exp (ikat) , 

k = ko + aki + ... 


where i fc 0 , h,.. ., 9 0 , #i, • • • are real constants, and h 0 ( t ) are periodic in time t. 

By writing the streamfunction and interfacial deflection in this form, the linear long wave 
stability of the unsteady flow can then be determined by calculating the first non-zero 9, 
which corresponds to exponential growth or decay of the disturbance to the basic flow. 

We now consider the system of governing equations (4a-b) and the boundary and interface 
conditions (5-11) to leading order, 0(1). Firstly the kinematic condition, (11) yields 

^ = -e 0 h 0 . (12) 

at 

Now since 9 0 is real and h 0 is periodic we have two choices: we may take h 0 = 0 so that $o 
is as yet unprescribed, or alternatively 6o = 0 and then, without loss of generality, ho 1 . 
We shall in fact assume the latter, since h 0 = 0 leads to damped modes only as shown later 
in Section 3.2. 

After making the transformation ipjo = <f>jo + Uj, for each region, j = 1,2, the leading 


order system becomes 


1 

O 

1 d 4 tp io 

dy 2 dt 

Re df ' 

1 

O 

<■0 

m 5 4 ^20 

dy 2 dt 

Re dy 4 


(y = M) = 1 + | [e int + C.C.] , 


(13a) 

(13b) 

(14a) 


6 


020 (y = 0,0 = 0, 


50io 

dy 

50 20 


(y = M) = 

(y = 0, t) = 


m 


am — m + 1 

1 


5y w ' am - m+ 1 
Here, c.c. denotes the complex conjugate. 


(3A 

2 

ts Ant 


— [L\e Pa - L 2 e /3 °) e‘ m + c.c. 


m2 


jKe + c.c. 


At the interface y = 1 we have 


5 3 0 


10 


010 = 020 , 

5010 5020 


-R, 


dy 

5 2 0io 
dy 2 
5 2 0 10 


= m 


dy 3 6 dy5f 

These equations admit solutions of the form: 


= m- 


5y ’ 
5 2 0 20 
dy 2 ’ 
5 3 0 20 


5y 3 


-B, 


5 2 0: 


20 


5yd<’ 


(14b) 

(14c) 

(14d) 


(15a) 

(15b) 

(15c) 

(15d) 


0io — Aio (y — a) 3 + Bio (y — a) 2 H — ~r + 1 

am — m + 1 

+ {[Cio sinh (/?y) + B 10 cosh (/?y) + E w y + B 10 ] e‘ n< + c.c.} , 

020 = A 20 y 3 + B 20 y 2 H — — — r 

am — m + 1 

+ | [C 20 sinh {f}ym~ 2) + B 20 cosh (fiym~ 2) + £ 20 y + F 20 ] e‘ nt + c.c.} . 


The twelve constants are found analytically by substitution into the above boundary and 
interface conditions, the resulting equations which determine A\q, A 2 o, ..., F 10, F 2 o are given 
in the appendix. 

At the next order, O (a), the kinematic condition (11) becomes 

dhi 

—77 — I- @1 + i0io (1,0 + ik 1 = 0 (16) 

Our solution strategy in obtaining the eigenvalues analytically, is based in part on the per- 
turbation eigenfunctions and interfacial amplitude being time periodic. In equation (16) 
above, 0io(l,<) is the sum of a real constant and a time periodic function, while k\, and 
#i are real. In order to obtain periodic solutions h\(t) the following constraints need to be 
satisfied, 

&i = —0io^ (1) = A10 (a — l) 3 — Bio (1 ~ a) 2 H — — : — 1, 

am — m + 1 
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( 17 ) 


hi = - J ixl>[ U o s) (l,t)dt, 

— -Q -1 [Cio sinh (/?) + Dio cosh (0) + E i0 + Fi 0 ] e ,m - c.c., 


where the superscripts ( 5 ) and {us) denote the steady and unsteady parts respectively. 

Before solving the 0(a) system, let us first consider the 0 (a 2 ) kinematic condition, 
which can be written as 

—7- — b @2 “I" fc&i ^1 “b 2 1 + i(j> 11 (1 5 £) = 0. 

at 


Using J to denote the imaginary part, we see that for to be periodic in t we require, 


0 2 = I 


{ 



Clearly then it is sufficient to solve for the steady part of the eigenfunction <j> u to 0(a) 
in order to determine 0 2 and hence the linear stability of the interface. We note that at 
this order, products of functions which have a time dependence given by exp (db'fH), yield 
additional steady terms in the governing equations and boundary/interface conditions. The 
eigenvalue 0 2 is therefore determined by both the steady and time oscillatory motion of the 
interface. The time independent perturbation equations are given by 


1 d 4 <f> {3) 


11 


Re dy4 
m d 4 <j$ 


(ki + Ui) 


rr \ d 2 <t> 10 d 2 U 1 


Re dy 4 

They admit a particular solution 


= 1 


( k 1 + U2 ) S . i 


dy 2 dy 2 

20 


dy 2 dy 2 


<f> 10 
<^20 


« 


d 2 u 2 l(s) 


r r [ a \( 1 ,7t \ ^io , ou i 

Xn = lRe k k k [^ + u ^~di~^ 0 ~di 


(*) 




X 21 = *m 


T7 \ 


du 2 l{s) 


dy 


dydydy , 
dydydy. 


(18a) 

(18b) 


(19a) 


(19b) 


so that the general solution of equations (18a-b) is 

<^n = (y - l) 3 + Bn (y - l) 2 + Cii (y - 1) + D u + Xn> 
(j> 2 i = A 2 i(y - l) 3 + B 2i (y - l) 2 + C 2 i(y - \) + D n + X 2 i- 


The eight constants are found by applying the following boundary and interface conditions: 


<f>n{y = a) = °> 


8 



4? 

a_M 

dy 

g4? 

dy 


(y = o) = o, 

(y = a) = 0, 


(y = o) = o, 




d<t> 


11 

,« 

11 


dy 

ay|? 

a s 2 

£4? 

a s 3 


— m- 


— m- 


+ Re 


$21 ) 

a4? 

dy 

.£4? 

<V 

a 3 4 ? 


(m-l) 


m 


Ai 


dU* 


lW 


dy J 


-/? e 


dy 3 

lk\d(f) 20 

dy 


ik\d<j>\o .jj d<f > io .j. dU\ 

+ iU 1-5 * 0 io 


|W 


+ if / 2 


dy 

d<t>20 


dy 


dy 

dU 2 


dy J 


dy 


(•) 


The above solutions and boundary conditions provide a system of eight inhomogeneous 
algebraic equations for the eight constants ^4n , ^ 21 , . . . , Z?n, D 21 and are solved explicitly 
here, (see appendix). 


3.2 Damped disturbances 

Before presenting the results of the calculation of the eigenvalue 0 2 , we show that the possi- 
bility h 0 = 0 yields damped waves. To achieve this we write the leading order streamfunction 
as 


<t>jo = Qj (y) exp (00 1) , 


where j = 1,2. 


This satisfies the equations 


which are solved subject to 


R e Q 0 

R e 0 0 


<PQ 1 

dy 2 

(PQ 2 

dy 2 


d 4 Q 1 
dy* ’ 

d 4 Q* 

m—r— 


<M o)-^(o) = 0 , 

0>(o) = ^W = 0 , 

<3i(i) = 9>(i), 

* 


( 20 ) 

( 21 ) 
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dQ '-(l) = fil(l), 


dy 


dy 


d 2 Q-i ^2 /lX 

^ (1) = 

_ ^r (1) ^-0) - m '^5- (1)_ Wo *‘ (1) - 

Multiplying (20) and (21) by ^ and ^ respectively, integrating between 0 and a , and 
imposing the above boundary and interface conditions yields 

' a {^y dy + k {^y dy\ 

R e 0 o = — m-Tr - ^ — — < 0. 


[/o (^) + fr (%) d y 

Hence all disturbances are damped for this case. 


4 Results 

The non-dimensional parameter A is the magnitude of the sinusoidal oscillations of the 
upper plate relative to the magnitude of the steady velocity. When A = 0, this problem 
reduces to that of steady plane Couette flow of two superposed fluids of differing viscosity, 
(as characterized by the ratio m). 

As an initial check on the validity of our analysis we firstly consider this steady case, which 
corresponds to the problem solved by Yih (1967). In order to make a direct comparison of the 
results we must first note that the length scale used by Yih (1967) is equivalent to (L - D) 
whereas in this work we use the lower fluid depth D . Yih also defines a parameter n = 
the depth ratio of the lower to upper regions, a brief calculation shows that n = (a — 1) 1 . In 
view of these notational differences the Reynolds number and streamwise wavenumber are not 
the same as their counterparts defined in this work, we must first make the transformations 
OLYih —3 ► na and Ry ^ nR e , we hence write J (m, a, A = 0) = R~ l (a — 1) 2 02- 

Calculation of 0 2 , and hence J, is in theory an analytical task, in practice however, the 
twenty constants obtained in the integration of the 0(1) and 0(a) systems, are lengthy 
expressions which are found most efficiently by use of a symbolic algebra package. For 
completeness, the simultaneous equations which determine these constants are given in the 
appendix. 

For a given viscosity ratio m, depth ratio a, frequency f I and magnitude A, the value of 02 
is calculated. When the real Floquet exponent 02 > 0 the disturbance to the basic flow will 
grow exponentially in time and the flow will be unstable, similarly negative 02 corresponds 
to a damped interfacial disturbance, and the flow is said to be linearly stable. 
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Figures 1(a) and 1(b) plot the values of J against m for a steady flow, (A = 0), when the 
lower fluid is more viscous than the upper layer. We have chosen equivalent depth ratios to 
those illustrated by Yih (1967), the results are identical. 

Figure 1(a) corresponds to a flow in which the upper layer is deeper, for a more viscous 
lower layer the growth rate is positive and the flow is always unstable. With a shallower 
upper layer, Figure 1(b) indicates that there is a region of stability, depending on the size of 
the viscosity ratio, although for sufficiently large m the flow will become unstable. 

For unequal depths of fluids the growth rate tends linearly to zero as m approaches unity, 
since this is a hidden mode for flows of equal viscosities. 

To obtain results for a less viscous lower region, we have verified the following transfor- 
mation given by Yih 


mJ (m, n, A = 0) = n 2 J f , A = (A, where n = (a — 1 ) 1 . 

\m n J 


( 22 ) 


The results for the steady flow when the more viscous fluid occupies the upper region can 
then be inferred from Figures 1(a) and 1(b). 

Before discussing the results for the full time-oscillatory problem, let us first investigate 
this property in more detail. Yih observes that the phase speed of the disturbed flow, 
(c(, using his notation), is equal in magnitude and opposite in sign when the viscosities 
and depths of the fluid layers are interchanged. This is because, the original flow may 
be recovered identically by making a Galilean transformation from the coordinates (x, y,t) 
to where ( = x - t, hence (£, £) -> A). The upper boundary 

now becomes stationary, (relative to the moving frame of reference), and the lower plate 
moves with constant velocity —1. Since gravity is neglected for this problem, inverting the 
geometry leaves the eigenvalue, and hence the stability characteristics, unchanged. In order 
to compare the growth-rates quantitatively however we must note that now the Reynolds 
number, R e = UoDp/pi, and lengthscale D must be rescaled on the new upper fluid, hence 
the multiplying factor n 2 /m, as given by equation (22) above. 

We now consider the problem when A ^ 0, the flow is dependent upon the time t and 
analogous conclusions from a change of reference frame do not hold and so for completeness 
we also need to consider arrangements when when the lower fluid is less viscous than the 
upper layer. Since the unsteady contribution is periodic, we need only consider A > 0. In 
fact it is clearly seen that A is a multiplicative factor in the basic flow, <^20 and also hi, 
and these terms contribute to $2 as products, thus —A gives identical results. 

When the lower fluid is more viscous, m > 1, the interface is destabilized by the intro- 
duction of an unsteady basic flow. Figure 2(a) shows that for a = 3.5 the already unstable 
steady mode, (broken line), is made more unstable upon increasing the magnitude of the 
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oscillations. Similarly, when the lower fluid is deeper, for example a = 1.4, so that in the 
absence of background time periodic modulations the flow is stable for moderate m at least, 
the interface is again destabilized as indicated by Figure 2(b). The value of m below which 
the flow is stable, decreases as A increases, for example when A = 0.2, the neutral distur- 
bance is obtained for a viscosity ratio m ~ 86.2319, whereas for A = 0.4, rn « 66.2506. This 
is shown more clearly by the neutral curve in Figure 3; if the magnitude of the oscillations 
increases beyond the critical value A 0.5695, the interface is unstable for all m > 1. 

When the lower fluid is less viscous, m < 1, the results are more significant. For all 
depth ratios we find that the time-dependent oscillations dramatically stabilize the interfacial 
disturbance. For a deeper upper layer, for example a = 2.25 as in Figure 4(a), the range of 
viscosity ratios m, for which the real growth rate is negative is increased for the unsteady 
background flow and as A is increased the flow is stabilized further. Figure 4(b), illustrates 
how the unstable steady mode corresponding to a shallow upper layer, can be completely 
stabilized provided A is made large enough. Such behavior is shown collectively in Figure 5 
which can be used to predict oscillation amplitudes which completely stabilize the flow for 
given depth ratios a. Figure 5 depicts the variation of neutral stability pairs m, A for three 
different depth ratios. In each case there is a global maximum oscillation amplitude above 
which the flow is linearly stable for all viscosity ratios. For example when the depth ratio is 
a = 1.4 we see that the flow is stable if A > 0.6412 (approximately). 

Finally we quantify the effect of varying the frequency of the oscillations, through the 
non-dimensional parameter U. Increasing the frequency reduces the effect of the unsteady 
terms and the stability of this fluid regime becomes comparable with that of steady plane 
Couette flow, (denoted by the broken line). These stability results can be understood by 
consideration of the unperturbed flow at high oscillation frequencies. In such a regime the 
flow separates into a Stokes layer in the vicinity of the oscillating wall, and away from 
this layer the flow is steady and corresponds to two phase Couette flow due to a boundary 
which moves with constant velocity. As long as the Stokes layer is thin compared to the 
thickness of the upper layer, therefore, the interfacial mode is expected to be insensitive 
to the wall modulations as indicated by our numerical results. The Stokes layer thickness, 
in nondimensional terms, is proportional to /?J 1/,2 fI -1/r2 , and for the highest values of f l 
used in Figures 6(a)- (b) the ratio between the Stokes layer thickness and the distance of 
the unperturbed interface from the wall is approximately 0.025 and 0.25 respectively for a 
unit Reynolds number, confirming the arguments given above. On the other hand, as fi is 
reduced the stabilizing or destabilizing tendencies characterized by A > 0, are emphasized, 
as shown by Figures 6(a) and 6(b). It can be concluded, therefore, that from a practical 
point of view, stabilization due to time harmonic modulations is likely to occur as long as 
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the modulation frequency is not too large. 

5 Conclusions 

We have considered the effect of the inclusion of time harmonic modulations in the unper- 
turbed velocity field of two phase Couette flow of different liquids. Our main aim is the 
demonstration of the stabilizing effect that such modulations can have on otherwise interfa- 
cially unstable flows, a finding that can have useful practical applications. To this end, we 
considered the unsteady partial differential stability system in the limit of long wavelength 
perturbations which can be solved exactly by use of Floquet and perturbation theory to 
yield analytical expressions for the first non-zero Floquet exponent which governs stability 
or instability. The main conclusion of this study is that inclusion of modulations with an 
amplitude larger than a certain parameter- dependent threshold, can completely stabilize 
flows which are unstable in the absence of modulations; for instance flows with the more 
viscous fluid occupying a thin layer and bounded by the moving wall can be completely 
stabilized for long waves, at least. At the same time, modulations can produce an adverse 
effect on the interfacial mode. Flow arrangements which are stable or unstable can become 
unstable or more unstable respectively. Finally, we emphasize that the conclusions of this 
work are valid for long wavelength perturbations which are useful in providing analytical 
solutions to a problem that needs to be addressed numerically in general. Further more, the 
physical conclusions reached for the long wavelength limit are expected to be indicative of 
the behavior of general wavelength time periodic perturbations. 


13 



Appendix 

For the solution of the 0(1) system (13a-b), we impose boundary and interface conditions 
given by equations (14a-d) and (15a-d). These yield the following twelve equations, which 
determine the constants Aio, A 2 q, ■■■, F 10 , F-io- 


A 10 — TnA'io 

3Aio (1 - a) - 3raA 20 + Bio — ™B 2 o 

A\o (1 — a) 3 — A 2 o + B\o (1 — a) — B 20 
• , m — 1 

3 A l0 (1 - a) 2 - 3 A m 4- 2 B 10 (1 ~ a) - 2 B 20 + — _ - + - 

E\o — E-20 

F\q — F 20 

E\o 4- E>2o 
A 

C\o sinh (0a) + O 10 cosh (/ 3a ) + E\oa + F\q — — 


0, 

0, 

0, 

0, 

0, 

0, 

0, 

0, 


[lae^ — L 2 eO^“) = ^Oio cosh (/3a) + ^Di 0 sinh(/3o) + £ 10 , 

pm 2 C 20 4- E 20 = P m 

C 10 sinh (0) + Dio cosh (0) = C 2 o sinh + D 20 cosh (0m~i) , 

Oi 0 cosh(/5) + D 10 sinh(/?) — [O 2 ocosh {0m~^ + O 20 sinh {0m 2 ) . 

At O (a) the solution of (18a-b), requires a further eight equations, to determine the constants 
An, A 2 i, Du, D 21 ' 


An ~ mA 2 1 = 0 

An (a - l) 3 + B11 (a - l) 2 + On (a - 1) + On + X11 (°) = 0, 
3An (a — l) 2 + 20n (a — 1) + On + (°) = 

— A 2 i + O 21 — 0 2 i + D 2 \ + X21 (0) = 0 
3A 2 i — 2 B 2 i + 0 2 i 4 — (0) = 0, 
On — 0 2 i 4 - X11 (1) — X21 (1) = 0, 


iT 7 1 (*) 


~ r ,5xii m 5X21 m (rra — 1) r dUi 


2B„-2mft l + ^(l)-m^l(l) = 0. 


Where Xn and X 21 are the particular solutions defined by equations (19a-b). 
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Fig. 2(a) Unsteady Effects: a=3.5, Q=1.0, m>1 
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Fig. 6(a) Frequency Effects: a=1.4, m>1 
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Fig. 6(b) Frequency Effects: a = 2.25, m<1 
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